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We study the physical implementation of the Photon Carnot engine (PCE) based on the cavity 
QED system [M. Scully et al, Science, 299, 862 (2003)]. Here, we analyze two decoherence mecha- 
nisms for the more practical systems of PCE, the dissipation of photon field and the pure dephasing 
of the input atoms. As a result we find that (I) the PCE can work well to some extent even in the 
existence of the cavity loss (photon dissipation); and (II) the short-time atomic dephasing, which 
can destroy the PCE, is a fatal problem to be overcome. 

PACS numbers: 03.65.Cw,05.70.Ca,42.50.Ar 



I. INTRODUCTION 

Recently many investigations have been carried out to 
explore various possibilities of building Carnot (or Otto, 
etc.) heat engines in some "quantum way" . It is expected 
that, by taking the advantages of quantum coherence, 
such quantum heat engines (QHE) using quantum matter 
as working substance can improve work extraction as well 
as the working efficiency in a thermodynamic cycle 0, 0, 
0, 3- Marian O. Scully and his collaborators proposed 
and studied a QHE based on a cavity QED system |j| 
0, Q, namely, a photon-Carnot engine (PCE) @. The 
working substance of their PCE is a lot of single-mode 
photons radiated from the partially coherent atoms. In 
their model, the walls of the cavity are assumed to be 
ideal, i.e., the cavity loss are disregarded. 

In practice, however, the walls can not perfectly reflect 
the photons (the cavity loss are not negligible), and the 
atomic dephasing of the input atoms are inevitable due 
to its coupling to the environment when passing through 
the cavity. To focus on the essence of the problem we 
only phenomcnologically consider the pure dephasing ef- 
fect |3 in this paper. A question then follows naturally: 
How does the photon dissipation and atomic dephasing 
influence the efficiency of the PCE? In this paper, by 
analyzing a more realistic cavity QED system, we re- 
vised the PCE model proposed by Scully and his col- 
laborators. We find the efficiency of the PCE decrease 
when the cavity quality Q becomes smaller (cavity loss 
becomes more strong); when the atomic dephasing hap- 
pens, though the atomic energy conserves, the quantum 
features of the PCE are demolished and then the QHE 
becomes a classical one. 

Our investigation is significant in two aspects. On the 
one hand, our results confirm the robustness of the PCE 
proposed in Refs. [a, Ifl 13j which can still work well 
even in the existence of not too strong cavity loss. On 
the other hand, our results demonstrate the quantum- 
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FIG. 1: (color online) The cavity QED model of our QHE. 
Three-level atoms, with quantum coherence in the two degen- 
erate ground states |<?i) and | c/2 ) , flow through the cavity and 
interact with the resonant photon field in the cavity. 



classical transition of the PCE due to quantum dephas- 
ing, which agrees well with our intuition, the efficiency 
of the PCE decrease due to atomic dephasing. These 
predictions can not only help us better the understand- 
ing of the basic concepts of thermodynamics, statistical 
mechanics, but also help us to optimize the system pa- 
rameters in future experiments of PCE. It is also of inter- 
est that the efficiency of the PCE in a Carnot cycle can 
be used to measure the quantum coherence of the input 
atoms and characterize the quantum-classical transition 
of the PCE. 



II. CAVITY QED MODEL OF QHE REVISED 

The PCE we consider here is similar to that proposed 
in Refs. 0, IS ( see the schematic illustration in Fig. 
I). In our PCE model, the ground states \gi) and 1 472) 
are accurately degenerate. The atom-photon coupling is 
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described by the Hamiltonian 0, [n| 



Hj = n\ \e) ( ■= )a + h.c. 

v2 



(1) 



Here, v is the level spacing between the excited state |e) 
and the ground states; a is the annihilation operator of 
the resonant photon field and A is the atom-field coupling 
constant. 

If there were no photon dissipation and atomic de- 
phasing, Hi would completely govern the pure quan- 
tum evolution. Let |m), m = 0, 1,2, ••• denote the 
Fock state of the photon field. In every invariant sub- 
space V m , which is spanned by the ordered basis vectors 
{|e) ® | m — 1) , \gi) ® |m) , \g2) ® \ m )}, the evolution op- 
erator U (r) = exp(— iHit) can be expressed as a quasi- 
diagonal matrix with the diagonal blocks 



U m (r) 



-^S m (T) —S m (^) 



V2 



"ml 2 / 
Cm(f ) 



(2) 



Here, C m (r) = cos (A^/mr) and «S m (r) = sin(AY / ror). 
In obtaining above explicit expressions for the matrix 
elements of U m (r), we have used the following technique. 
By writing 



\G) 



(.91 1 + (.92 1 
V2 



we find 



Hi 



fiA|e) (G|a + ft.c 



(3) 



is, in fact, the Hamiltonian for the Jaynes-Cummings 
model at resonance, and then we can exactly solve the 
problem of time evolution in each subspace V m following 
the method in Refs. [H [n). 

We denote the initial density matrix of the photon 
field and the atom ensemble as p L (ti) and p A (ti). The 
reduced density matrix p L (t) of the radiation field will 
evolve according to the "super-operator" M (r) defined 
by 

M (r) p L (ti) = Tr A [U (r) p L (U) ® p A (U) (r)], (4) 

where Tr A means tracing over the atomic variable. The 
atoms pass through the cavity at the rate r. Then we 
can write the known master equation at zero temperature 
[H[lll as 



-p L (t)wr[M(r) 



l]p L (t) + Lp L (t) . 



(5) 



Here, we have made the approximation In [M (r)] w 
M (r) — 1 for a short time r (the reliability of this ap- 
proximation can be seen from later parameter estimation 
r ~ 10~ 10 s), and the cavity loss term is defined as 

Lp L (t) = (^)[2a Pi (i)at - {o)ap L (t) + h.c.)] (6) 



with the cavity quality factor Q. 

Inspired by the concept of "phaseonium" in Ref. , we 
prepare the atoms initially in a partially coherent state 



PA (0)= Pe \e) (e\ + \g) (g\. 



(7) 



Here, \g) (g\ contains the superposition of the ground 
states 

Iff) (.91 = |ci| 2 k)(ffi| + |c 2 | 2 |5 2 )(ff 2 | (8) 
+cic* 2 \gi) (52 1 + c\c 2 \g 2 ) (gi\ ; 

p e is the probability distribution in the excited state |e). 
p A (0) has been normalized to unit, i.e., 



Pe + \Cl\ 



\C2\ =1- 



(9) 



It is well known that because the atomic dephasing is 
much more rapid [ill IT2I Il3| than the quantum dissipa- 
tion of the atoms, quantum dephasing is the dominating 
decoherence mechanism in short time. Hence, we only 
phenomenologically consider the atomic dephasing when 
the atoms pass through the cavity. We replace the pure 
state \g) (g\ in p A (0) with the mixture 

Pg = E ' Cfc | 2 ^ ^ + ^ ClC 2 \9i) (521 + h.c, (10) 

fe=l,2 

where ^ is the so-called decoherence factor satisfying |^| < 
1 0. The complete coherence is characterized by £ = 
1, while the complete decoherence is by £ = 0. With 
above considerations, we depict the atomic dephasing by 
changing the initial state from p A (0) to 



Pd = Pe 



(11) 



III. PHOTON- "WORKING FLUID" AT 
STATIONARY STATE DESCRIBED BY AN 
EFFECTIVE TEMPERATURE 

Next we consider a similar cycle as that in Ref. 
0. During the isothermal expansion process the in- 
put atoms, prepared in a state slightly deviating from 
the thermal equilibrium with some quantum coherence, 
serves as a high temperature energy source, while during 
the isothermal compression process the input atoms, pre- 
pared in thermal equilibrium state, serves as the entropy 
sink. For a short interaction time r, we have (the relia- 
bility of this approximation can also be seen from later 
parameters estimation y/mXr ~ 1CP 1 ) 



C m (r) ~ 1 — -m\ 2 T 2 , 



(12) 



S m (r) 



mXr. 



Then we obtain the equation of motion for the average 
photon number (n (t)} 

I („ (t)) = M [(2 Pe - 9) (n (i)) + 2 Pe ] -Un (t)) , (13) 
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where 



2^2 



(14) 



\c 2 \ 2 + 2Re(^c 1 c* 2 ), 



and Re(£,cic 2 ) is the real part of £,c\c 2 . 

In the thermal equilibrium state, the atomic probabil- 
ity distributions p e , |ci| 2 and | C2 1 2 satisfy 



C2 



Til/ 

ex P (-^), 



(15) 



where A: is the Boltzmann constant and T is the tem- 
perature of the thermalized atoms. Since the relaxation 
time of the radiation field is very short, in the following 
analysis, we will use the equilibrium state solution 



{n E ) = 



1 + CCO 



(16) 



of Eq. Ijl3|l to replace the time-dependent average photon 
number (n (t)). Here, n is the average photon number 



2 Pe 



\c 2 \ 2 -2p e 



(17) 



in the absence of atomic coherence and cavity loss, and 



C(T) = -[Re(£ Cl ^) , 



(18) 



is a temperature dependent parameter concerning the 
cavity loss as well as the atomic dephasing. 

We imagine the radiation field also obeys a virtual Bosc 
distribution 



1 



{nE) = 



exp[hu/(kT')] - 1 



(19) 



with an effective temperature T'. In high temperature 
limit, we approximately have (for the microwave cavity 
QED system and circuit QED, this approximation is re- 
liable when the effective temperature T' is at room tem- 
perature or higher) 



{n E 



kT' 



kT 

hv 



(20) 



Therefore T' can be approximately determined as |j| 



T' = 



T 



i+ccn 



(21) 



It can be seen that the effective temperature T 1 being dif- 
ferent from T is due to the atomic coherence as well as the 
cavity loss. Obviously, when Q — > oo and Re^cicJ) = 0, 
the effective temperature T 1 approaches T. I.e., the ef- 
fective temperature becomes equal to the temperature of 
the input atoms when the atomic coherence cancels and 
the cavity is perfect. 
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FIG. 2: (color online) Temperature-entropy diagram for the 
photon-Carnot cycle. Here we consider neither the cavity loss 
nor the atomic dephasing. i.e., £ = 1, Q — > oo. In addition, 
we consider the special case of the phase angle Arg(ciC2) = n. 
As a result the effective temperature T' h is higher than Th by 
the amount AT h = n h T h \c\c\*\IP^. 



IV. THERMODYNAMIC CYCLE OF THE QHE 
WITH QUANTUM MATTER AS THE WORKING 
SUBSTANCE 

The Carnot cycle of our QHE consists of two isother- 
mal and two adiabatic processes @ (see Fig. 2). We 
use the subscripts (or superscript) h and I to indicate 
the isothermal expansion and the isothermal compres- 
sion processes hereafter. E.g., and ni represent n in 
the isothermal expansion and the isothermal compression 
processes respectively, i.e., 



2PI 



J 12 



„i 12 



2p% 



hj. 



(22) 



During the isothermal expansion process, from a thermal 
state 1 of the photon field to another 2, the input three- 
level atoms are prepared with quantum coherence of the 
ground states p h A (0) = p A (0). But during the isothermal 
compression process from 3 to 4, the input atoms are 
prepared in a thermalized state, i.e., 

Pa (0) = Pi \e) (e| + |4| 2 \ 91 ) (g x \ + |4| 2 \g 2 ) (g 2 \ . (23) 

^From above facts and Eq. (|18fl we know, in the isother- 
mal expansion process from 1 to 2 



C h (T h ) = ^[Re^c14*) 

Pe 



— ] 



but in the isothermal compression process 
r (T\ - vm 

Cl{Tl) -2-M- 

We apply the entropy expression 

hvi (n l E ) 



(24) 



(25) 



S t = k\n((n E ) + l) + 



T! 



(i = h, I) (26) 



of the radiation field to calculate the heat transfer. In 
a Carnot circle, the work done by the radiation field is 
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AW = Qi 



and 



Qout- Here 
Q in = T' h [S h {2)-S h {l% 

Q ou t = Tl[S l (3)-S l (4)} 



(27) 



(28) 



are respectively the heat absorbed into the cavity dur- 
ing the isothermal expansion process from 1 to 2 and 
the heat released out of the cavity during the isothermal 
compression process from 3 to 4. 

Similar to the cycle in Ref . , the frequency of the 
radiation field, i.e., the resonant mode of the cavity, is 
assumed to change slightly form 1 to 2, i.e., 



HI) -^(2)] 
u{l) 



< 1. 



Namely, we can make the approximations 



1/(1) 
1/(3) 



Vh 



"(2), 
<i/(4). 



(29) 



(30) 



In the adiabatic process, the average photon number does 
not change, i.e., 



K (2)) 
(n h E (4)> 



(n l E (3)) 
{n l E (1)) 



(31) 



Form Eq. ("2""J), it follows that 

1/(1) 1/(4) i/(2) 



7" 



IT 



T" 



"(3) 



(32) 



From these observations and Eq. (|26|) we find that 

S h (2) - S h (1) = Si (3) - S, (4) . 

Therefore, in the high temperature limit, the PCE effi- 
ciency rj = (Q in — Qout) IQin can be explicitly expressed 
as V = 1 - T[IT h I or 



*7 



1 



1 + Ci(Ti) 



(33) 



Based on above results we are now able to analyze the 
effects of the two dccohercncc mechanisms separately. 

Firstly, to focus on the cavity loss (photon dissipation) , 
we consider the ideal case with no atomic dephasing, £ = 
1. In the ideal case the cavity loss is negligible, i.e.,Q — > 
00, the efficiency (|5*3*|) is reduced to 



V 



1-1 



Pe T h 



(34) 



which agrees with the result in Ref. [Bj]- It seems that, in 
principle, the PCE can extract positive work from a sin- 
gle heat bath if we control the phase angle 8 =Arg(c^C2*) 
properly, e.g., 8 = n. This shows the advantage of 
the "quantum fuel" -we can extract more work from the 



"quantum fuel" than from the classical fuel. However, 
in the case the cavity is "extreme bad", i.e., the cavity 
quality factor is vanishingly small Q — > 0, the efficiency 
decrease to zero (this can be verified easily) 



n 



1 



Vhn h P e Ti 
viniph T h 



0. 



(35) 



Namely, the PCE is destroyed by the strong loss of the 
cavity. This can also be understood intuitively from the 
Eq. i|2Tjl. When the quality factor of the cavity becomes 
so small that few photons can stay stably in the cavity. 
Accordingly, both the two effective temperatures T' h and 
T[ becomes vanishingly small, and thus no work can be 
done by the "working substance". As a result, the effi- 
ciency of the PCE decrease to zero. Mathematically, it 
can be verified that r\ is a monotonically increasing func- 
tion of Q, and the efficiency of the PCE decrease to zero 
when Q becomes vanishingly small. 

Secondly we consider the pure atomic dephasing. From 
Refs. 0, 0, we know that the dephasing time is 
much shorter than the atom and cavity lifetimes. After 
the atom interacting with the environment for a short 
time t, the term Re^c^c^*) concerning atomic coherence 
becomes vanishingly small. We can properly assume £ 
decrease to zero in the short time r. Then the efficiency 
i|53| of the PCE becomes 



7J = 1- 



1 + v h n h /{2^p h e Q) 
1 + vim/(2np l e Q) 



(36) 



In principle, if the cavity is ideal, Q — > 00, we regain the 
maximum classical Carnot efficiency r\ = 1 — T\ /Th from 
the efficiency (|36|l of the PCE. It turns out that, without 
cavity loss, the complete dephasing of the atoms makes 
the PCE become an ideal (reversible) classical heat en- 
gine. This demonstrate the quantum-classical transition 
due to the atomic dephasing. Similarly, in the case of 
"extremely bad" cavity, Q — > 0, the efficiency (|36|) de- 
crease to zero (15*51) . 

Finally we analyze the positive- work condition of this 
PCE [T4I "T5L IT""(. under which positive- work can be ex- 
tracted. From Eq. Q33[l we know that the positive work 
condition AW = Qi n — Q ou t > of this QHE is 



lh > i+m) h 



(37) 



where [l+CftP/i)]/[l+Ci(-^)] can be either less or greater 
than unit. It is counterintuitive when [1 + £ h (Th)]/[l + 
Ci(Ti)] is smaller than unit, and the same novel result 
occurs in Ref. {5J. These novel results originate from 
the fact that the atoms are not prepared in the thermal 
equilibrium state in the isothermal expansion process. In 
other words, the initial state with partial quantum coher- 
ence is out of thermal equilibrium, and thus the "temper- 
ature" Th of the input atoms in the isothermal expansion 
process should not be regarded as a real thermodynamic 
temperature essentially. We also would like to emphasize 
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FIG. 3: (color online) Temperature-entropy diagram for the 
photon-Carnot cycle. The atomic coherence are demolished 
due to the atomic dephasing, i.e., f = 0. The cavity loss 
makes the effective cavity temperature decrease. Here, ATh = 
n h T h v h / (2/iPjQ), and AT, = mTiVi/ (2/iP^Q). 

that the second law of thermodynamics is not violated 
when 

for it would take energy from an external source to pre- 
pare the atomic coherence 0. In an overall considera- 
tion, the extra energy cost for preparing the atomic co- 
herence would prevent the second law from being vio- 
lated. 



Vhj (2/xQ) in Eq. I|24|) has the order of magnitude greater 
than 10 _1 , but the order of the other part Re(fcyc2*) in 
Eq. i|24|) is less than 10 _1 due to the atomic dephasing 
£ < 1. Thus the conclusion for optical cavity system 
remains valid for microwave cavity system. 

Thirdly, we consider the circuit QED based on super- 
conducting Josephson junction systems [ItI Hsj . Here, 
the coupling between the charge qubit [a Cooper pair 
box (CPB)] and a quantum transmission line is in the 
similar way as the photon-atom coupling in cavity QED 
system. The passage of atoms through the cavity can 
be simulated by the periodic switch -off and -on of the 
on chip- interaction |l9l l2(i| . The experimental parame- 
ters can also be found in Ref. 0: the resonance fre- 
quency Vh r~j 10 10 Hz; the atom-held coupling con- 
stant A ~ 10 s Hz and the quality factor is less than 
10 4 . After similar analysis, we hnd the value of the term 
Vfi/ (2/iQ) in Eq. I|24|) has the order of magnitude greater 
than 10°, but the order of magnitude of the other part 
Re(£ciC2*) in Eq. J23J is less than 10 _1 due to the atomic 
dephasing £ < 1. Therefore, under current experimental 
capability, the Circuit QED system can not also be used 
to implement the QHE, either. 

We therefore conclude, based on the present experi- 
mental accessibility, we will have to improve the quality 
factor of cavity Q to a much higher level before we can 
implement the PCE as a practical QHE. 



VI. CONCLUSIONS 



V. EXPERIMENTAL FEASIBILITY 

Finally, we would like to estimate the efficiency rj ac- 
cording to a set of experimentally accessible parameters. 
For different material and cavity parameters [ITI the dis- 
cussions are listed below. 

Firstly, for the optical cavity QED system ^?|? the 
resonance frequency (vi) ~ 10 Hz, the atom-field 
coupling constant A ~ 10 8 Hz and the quality factor 
is less than 10 8 . We take a set of reasonable values 
H @: n h (n t ) ~ 10 2 , Xt^/E ~ 1CT 1 . Thus we get or- 
der of magnitude of the interaction time r ~ 1CP 10 s. 
To require only one atom in the cavity once, we have 
r < (1/ T ) ~ 10 10 /s. Based on above estimation, we get 
fx = rX 2 r 2 /2 < 10 6 /s. Hence, the term v h j {2fiQ) in Eq. 
(|24H has the order of magnitude greater than 10°, but the 
order of magnitude of the other part Re(£cjc2*) in Eq. 
(|24|l is less than 1CP 1 due to the atomic dephasing £ < 1. 
It can be seen, for the optical QED system, when consid- 
ering the practical experimental parameters, we can well 
approximately disregard the term Re(£cYC2*) concerning 
quantum coherence in calculating the efficiency. 

Secondly, for the microwave cavity QED system [T^| . 
the resonance frequency Vh (vi) ~ 10 10 Hz, the atom- 
field coupling constant A ~ 10 4 Hz and the quality factor 
is less than 10 9 . After similar analysis, we find the term 



In summary, we revised the PCE proposed in Ref. 0, 
and we found that the PCE efficiency decrease monoton- 
ically with Q due to the cavity loss. In the ideal case, 
Q — > oo, we regained the result of Ref. |5(, where the 
improvement of the working efficiency is due to the non- 
equilibrium state preparation. Thus the obtained effi- 
ciency, which is beyond the classical Carnot efficiency, 
does not imply the violation of the second law of ther- 
modynamics. This observation has been made in Refs. 
[l4l Il5l | for different physical system. We also phe- 
nomenologically considered atomic dephasing and found 
that a short time dephasing may make the PCE become 
a classical one. From these heuristic discussions, we con- 
clude that both the photon dissipation and atomic de- 
phasing can diminish the efficiency of the QHE. Con- 
sidering the continue improvement of the cavity quality 
factor Qi we believe the crucial issue in future experi- 
ments will be to control the atomic coherence. We also 
would like to point out that the dissipation mechanism 
of atoms due to its coupling with the environment, e.g., 
the vacuum modes, are not considered microscopically 
in this paper. Detailed investigation of the PCE with 
atomic dissipation will be presented in our forthcoming 
paper. 
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